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Quantised vortices are remarkable manifestations on a macroscopic scale of the coherent nature of 
quantum fluids, and the study of their properties is of fundamental importance for the understanding 
of this peculiar state of matter. Cavity-polaritons, due to their double light-matter nature, offer 
a unique controllable environment to investigate these properties. In this work we theoretically 
investigate the possibility to deterministically achieve the annihilation of a vortex with an antivortex 
through the increase of the polariton density in the region surrounding the vortices. Moreover we 
demonstrate that by means of this mechanism an array of vortex-antivortex pairs can be completely 
washed out. 

PACS numbers: 03.75.Lm, 42.65.Hw, 71.36.+C 


I. INTRODUCTION 

Topological excitations such as quantised vortices, 
characterised by a phase winding from 0 to 27rm (with 
m an integer number) around a vortex core, have 
been extensively studied in several systems such as 
non-linear optical system^, superconductor^, superfiuid 
vertical-cavity surfac e-emitting laserJ^, and more 
recently in cold atomJ^J^. Einally, in recent years, the 
study of vortices and vortex lattices has attracted much 
attention also in the field of coherent cavity-polariton flu¬ 
ids. Eirst because, being intrinsically out-of-equilibrium, 
they constitute a novel system to study Bose-Einstein 
condensation phenomena, and secondly because, due to 
their light-matter nature, polaritonic systems are fully 
controllable by optical techniques and therefore allow 
very detailed studies of quantum turbulence. 

In particular, in the context of cavity-polariton sys¬ 
tems, i t has been shown that stable vortices and half- 
vortice^^^^ as well as single vortex-antivortex (V-AV) 
pair^nmUMlI be generated. The formation of lat¬ 
tices of vortices and of vortex-antivortex pairs has also 
been theoretically and experimentally studied in sev¬ 
eral different configurations: in the optical paramet¬ 
rical oscillator co nfignrati orP^, in non-resonantly gen¬ 
erated condensateand in the case of patterns 
induced by metallic deposition on the surface of the 
cavit^ff^. However, the mechanisms lying beneath vortex- 
antivortex annihilation and vortex-vortex interaction are 
still not fully understood. This is related to the high 
degree of control needed to study such interactions. 

Eor example, in the fir st exper imental observations of 
vortex-antivortex lattices^^^^^ the formation and the 
properties of the array were only partially controllable 
due to either the presence of a strong exciton reservoir. 


that influences the position of the formed vortex array 
and its disappearance, or due to the fact that the formed 
array depends on the structure of the metallic depositions 
over the cavity surface. To achieve a higher degree of 
control resonant pumping schemes have been proposecP^. 
By using masks in the pumping beam, the formation of 
vortex-antivortex arrays with controllable shape and vor¬ 
tex distribution was achieved^ and the evolution of the 
arrays was studied in correlation with the local onset of 
the superfluid regime. However, in this study the case of 
high polariton densities was not experimentally achiev¬ 
able since the masks used to generate the vortex array 
were blocking most of the laser power, and therefore was 
not theoretically investigated. 

In this work we theoretically study the annihilation of 
a vortex with an antivortex when injecting polaritons res¬ 
onantly, in a broad range of polariton densities. In our 
model we assume an excitation with four coherent laser 
beams resonant with the lower polariton branch similarly 
t(P3an(P^ where, however, polaritons were injected non- 
resonantly. The four pump spots are supposed to have 
the same energy and k-vector modulus and to generate 
polaritons propagating toward the centre of a common 
area. The advantage of our model is that in the reso¬ 
nant configuration it is well known that there is no ex¬ 
citon reservoir and that, since the entire pump intensity 
is used to inject the coherent fluid, high polariton densi¬ 
ties can be achieved. Moreover, since the pumps set the 
momentum and the density of the injected fluid we can 
directly correlate the vortex-antivortex annihilation and 
the washing out of an array of vortices with the increase 
of the regions where the fluid is subsonic and with the 
change in the polariton flow that comes with it. 

The manuscript is structured as follows: In Sec. [ll|we 
introduce the theoretical model used to simulate the sys- 
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tern, describe more in details the setup we have in mind 
and introduce the generalized loeal speed of sound that 
will be needed for the analysis of the results. In Sec. 


Ill we demonstrate the annihilation of V-AV pairs and 


correlate this annihilation with the change in the polari- 
ton flow caused by the widening of the regions where the 
fluid is subsonic. Finally in In Sec. |IV| we draw some 
conclusions and give prospects for future developments 
of the work. 


II. MODEL 

A standard way to model the dynamics of resonantly- 
driven polaritons in a planar microcavity is to use a 
Gross-Pitaevskii (GP) equatioiP^ for coupled cavity and 
exciton fields and generalized to include the 
effects of the resonant pumping and decay {h= 1): 
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where the single particle polariton Hamiltonian Hq is 
given by 

^ \ ^i?/2 LUc {—i^) — ii^cl‘^) ' 


and 


ujc — ^c(O) 


2mc 


is the cavity dispersion, with the photon mass me = 
5 X 10“^mo and mo the bare electron mass. For our 
simulations we assumed a flat exciton dispersion relation 
co’x(k) = Cc;x(0), set the exciton-photon detuning to zero 
^ex-ph = ^x(O) — ooc{^) = 0 and set this energy value 
as the reference of the zero energy. The parameters Qr, 
hzx and hzc are the Rabi frequency and the excitonic and 
photonic decay rates respectively and have been given 
values close to the usual experimental ones: Qr = 5.1 
meV, Kx = 0.05 meV, and i^c = 0.08 me\^. In this 
model polaritons are injected into the cavity by four co¬ 
herent and monochromatic laser fields with pump inten¬ 
sity fp and Gaussian spatial profiles with ap of 20 jam: 

Where kp. are the 
four wave vectors of the four pumps that we fix to have 
the same modulus |kp. | = |k|. In order to ensure that 
the phase of the polariton fluid is not imposed by the 
laser pumps in the central region of the system we set to 
zero the pump intensity outside of a apin = 9/im radius 
circle. The exciton-exciton interaction strength gx is set 
to one by rescaling both the cavity and excitonic fields 
and the pump intensities. The numerical solution of the 
GP equation is obtained over a two-dimensional grid (of 
512 X 512 points) in a box with sides of 150 x 150 
using a fifth-order adaptive-step Runge-Kutta algorithm. 


All the analysed quantities are taken when the system has 
reached a steady state condition after a transient period 
of 200 ps. 

To understand the role of the subsonic character of 
the fluid in the annihilation of a vortex-antivortex pair 
we define the local fluid velocity = fi|k(x)|/mLP, 

where mRp is the lower polariton mass and k(x) is 
the locally evaluated derivative of the phase at the 
point X. Moreover we define the quantity Cs(x) = 
yfigLp|4'Lp(x)|2/TOLP, where |5 'lp(x)| 2 and qlp are 
the local polariton density and the coupling constant. 
Since in the local density approximation Cs(x) corre¬ 
sponds to the speed of sound defined in the case of high 
densitieJ^^^^, we can take it as definition of a generalized 
loeal speed of sound valid also for low polariton densities 
and we define a generalized Mach number: 


M(x) = = h\k{x)\/mLP 

Cs (x) \/fiflLp|^'Lp(x)p/TOLp 

This will allow to establish a direct correlation between 
the subsonic or supersonic character of the fluid and the 
annihilation of a V-AV pair. Since the subsonic charac¬ 
ter of the fluid is induced by polariton-polariton interac¬ 
tion, this corresponds to studying the role of polariton- 
polariton interactions in the annihilation of the pairs 
and of the vortex-array. Although other techniques, like 
the study of the vortex-antivortex correlation functiorP^ 
or Reynolds-averaged Navier-Stokes equationP^, can be 
used to address this problem, we chose to focus on the 
study of the Mach number since it since it allows a simple 
and clear physical understanding of the V-AV merging 
process. 


III. RESULTS 

The mechanism lying beneath the annihilation of V- 
AV pairs and the role of polariton-polariton interactions 
in this annihilation can be better highlighted by studying 
the system behaviour as a function of the pump inten¬ 
sity and therefore of the polariton density. We start by 
studying the system in the two limiting cases of very low 
and very high pump intensities. 

In the low intensity case the polariton density lies on 
the lower branch of the bistability curve everywhere in 
space and the system behaviour is purely linear. In this 
regime (fig. Eh -G), the formation of an array of vortices 
and antivortices is observed as iiinmUIl. Here, like iiP, 
the shape and size of the unit cell only depend on the 
geometry of the pumping configuration and on the angle 
of incidence of the laser beams. Since we use four pumps 
with |k| = 0.7 the formed array has square unit cells 
with unit cell size of approximately 9/im (fig. [T- G).This 
interference pattern generates an array of vortices (with 
a clockwise phase winding from —tt to tt) and antivortices 
(with an anti-clockwise phase winding from — tt to tt) that 
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FIG. 1: Numerical Real space emission intensity, phase, 

and vortex distribution in the low and high density regimes. 
The four polariton ensembles are generated by pumps at 
kpump = 0.7 /im“^ in different directions (yellow arrows) with 
(j^pump = —2.25 meV. A: Real space image in the low density 
regime with pump intensity fp = 0.07fth, a square interfer¬ 
ence pattern with a unit cell of about 9/im is formed. D: real 
space image in the high density regime fp = 3.33fth- B and 
E: Phases corresponding to A and D showing an array of vor¬ 
tices and antivortices (B) and no phase modulations (E). C 
and E: Vortices (red) and anti-vortices (blue) distribution in 
real space corresponding to phase diagrams B to E. 


is therefore due to purely linear mechanism (fig. wC and 
inset for the definition of vortices and antivortic^. 

In the opposite limit the pump intensity is strong and 
polariton-polariton interactions play a dominant role. 
This regime is characterised by a polariton density every¬ 
where lying on the upper branch of the bistability curve. 
As it can be seen in fig. EP -F, in this regime the array of 
vortices and antivortices has completely disappeared due 
to the renormalisation of the lower-polariton branch. 

From this, we can deduce that in the transition region 
between the low and the high density regimes vortices 
and antivortices either are expelled from the fluid or an¬ 
nihilate each-other. In order to investigate this transition 
we vary the pump intensity around the threshold value 
ifth) a.t which V-AV pairs disappear. Note that in this 
non-homogeneous system the density if different in the 


regions within and outside the pumping spot, resulting 
in four threshold intensities: two for increasing and two 
for decreasing pump intensities. When the intensity of 
the four pumps increases the regions directly pumped 
by the lasers jump from the lower to the upper branch 
of the bistability curve (first threshold). At this point, 
since polaritons have a finite lifetime, the central region 
between the four laser spots is still in the lower part of 
the bistability curve. As the pumps intensity is further 
increased also the central region eventually jumps from 
the lower to the upper branch (second threshold). This 
threshold corresponds to the intensity fth at which V-AV 
pairs disappear. Similarly one can observe two threshold 
for decreasing pump intensity. 

Figures A-C(D-F) represent the polariton distribu¬ 
tion (phas^ for increasing pump intensity from just be¬ 
low to just above fth- In fig. [2|V {fp = OMfth) some 
remaining of the interference pattern of fig. [T}\ is still 
visible together with four dark segments surrounding the 
centre of the image that correspond to four V-AV pairs. 
In this plot of the intensity distribution each V-AV pair 
looks like a straight dark segment rather than like two 
separated vortices because the core of the vortex is ex¬ 
tremely close to the core of the antivortex. The fact 
that these four dark segments correspond to V-AV pairs 
is confirmed by the phase distribution of fig. IP where 
the two phase rotations of the vortex and of the antivor¬ 
tex are visible in correspondence of each dark segment 
of fig. |2}\. For the sake of clarity the phase distribution 
of the V-AV pair delimited by the red square in fig. IP 
is also reported enlarged in fig. IP- When the pump in¬ 
tensity is increased (fig. (fp = O.Sfth) and fig. IP 
(fp = 1.0fth)) the vortex and the antivortex cores get 
closer (i.e. the dark segments become shorter) until the 
four V-AV pairs disappear. For even higher pump inten¬ 
sities the density distribution becomes homogeneous as 
in fig. IP- Again this moving closer and merging of the 
V-AV is confirmed by the corresponding phase distribu¬ 
tions (fig. IP -F), and by the corresponding zoom of the 
regions delimited by the red squares (fig. Ip-i). 

This detailed analysis of the density and phase distri¬ 
butions and of the ”on site” annihilation of V-AV pairs 
shows that the disappearance of the vortex array is not 
due to the expulsion of the vortices from the fluid. There¬ 
fore, the mechanism lying beneath the V-AV annihilation 
cannot be ascribed to a simple renormalisation of the 
lower polariton branch. First, the effect of the renor¬ 
malisation is to decrease the wave vector of the injected 
polaritons and correspondingly to increase the size of the 
interference pattern of the array therefore leading to the 
expulsion of vortices from the fluid rather than the ob¬ 
served merging. Second, it can be seen in fig. IP that 
when V-AV annihilation takes place, part of the interfer¬ 
ence pattern is still visible. Third, vortices and antivor¬ 
tices can exist in a fluid at rest, so the renormalization 
of the lower polariton branch can not justify per se the 
disappearence of the array. 

Moreover, since all the plots show the system steady 
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FIG. 2: A-C: Numerical Real-space emission intensity for 
three different pump intensities. D-F: real space phase- 
diagrams corresponding to the pumping condition of A-C. 
G-I: Enlargement of the red squared area in D-F. L-N: Mach 
charts corresponding to panels A-C, where the fluid is super¬ 
sonic the Mach chart show the real-space emission intensity 
at that point, when the fluid is subsonic the Mach chart is 
blue. A V-AV pair is clearly visible in G, V and AV have 
almost merged in H, and no vortex pair is present any more 
in I. The pump parameters are the same as in fig. [^and the 
pump intensities (from left to right) are fp = 0.66/t/i, 0.8fth, 
l.Ofth. 


state for a given pump intensity, the fact that vortices 
coexist with antivortices means that their position is the 
result of a tradeoff between vortex-vortex and vortex- 
antivortex interactions, direction of the polariton flow, 
and relative phase of the four pumps. In other words, 
this means that vortices and antivortices are not com¬ 
pletely free to move in the system and therefore V-AV 
annihilation cannot be simply ascribed to V-AV attrac¬ 
tion. In fact, if vortices and antivortices are completely 
free to move, due to their mutual attraction they must al¬ 
ways annihilate, and the only possible steady state must 
be completely free of vortex dislocations independently 
from the intensity of the laser pumps. Finally, it is worth 
noting that due to the choice of the continuous-wave reso¬ 
nant pumping setup no exciton-reservoir is present in the 
system and all polaritons have the same energy therefore 
no trapping mechanism can be advocated to explain this 
annihilation as iiP. 

To have a physical insight of the mechanism lying be¬ 
neath V-AV annihilation it is useful to study the Mach- 



FIG. 3: A-G: Numerical real-space emission intensity (in lin¬ 
ear greyscale) . D-F: Real space phase diagram (color scale) 
corresponding to A-G. The red squares indicate the regions 
where vortex and anti-vortex pairs have been annihilated. 
The plots correspond to pumps with increasing in plane mo¬ 
menta from left to right: kpump = 0.5, 0.7 and 0.9 /im“^. In 
each case the pump intensity is chosen in order to drive the 
system at the threshold at which V-AV annihilation occurs 
{fp{B) = 4.33 fp{A)Jp{C) = 8.66fp{A)). 

number charts in fig. |^-N where blue regions correspond 
to a subsonic character of the fluid and grey regions cor¬ 
respond to a supersonic character. Clearly as the pump 
intensity is increased (from left to right) the polariton 
density increases and the regions where the fluid is sub¬ 
sonic becomes wider and wider. This spreading of the 
subsonic regions rearranges the polariton flow therefore 
changing the steady state position of vortices and an¬ 
tivortices until, when the pump threshold fth is reached, 
vortices and antivortices annihilate. This mechanism can 
also be understood observing that vortices cannot en¬ 
ter subsonic regions since they are intrinsically related 
to regions where the fluid is supersonic, because at the 
centre of their core the fluid density vanishes while the 
fluid velocity does not and because a subsonic fluid tends 
to avoid strong phase modulations^. Since in our sys¬ 
tem the subsonic regions surround the V-AV distribution, 
when the pump intensity changes the polariton flows re¬ 
arrange letting vortices overlap and annihilate with an- 
tivortices. This is different from the case of atomic Bose- 
Einstein condensates where a single vortex cannot be de¬ 
stroyed by the superfluid character of a quantum fluid 
since the total angular momentum of the system must 
be conserved. Here, polariton-polariton interactions can 
affect the fluid distribution and make a vortex overlap 
with an antivortex, so that V-AV pairs can be completely 
suppressed in the fluid. 

An additional proof of the correlation between the sub¬ 
sonic character of the fluid and the disappearance of V- 
AV pairs is given by the fact that higher sound velocities 
(i.e. higher polariton densities) are needed to destroy 
the array when polaritons have higher velocities. Fig¬ 
ures A-C show the real space distribution of the pho¬ 
tonic field for three different increasing velocities of the 
injected polaritons {kpump = 0.5, 0.7 and 0.9 /im“^) and a 
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pump intensity corresponding to the threshold {fp = fth) 
where V-AV pairs disappear. We find that the polariton 
density, and therefore the sound velocity, at which V-AV 
pairs annihilate is higher when the velocity of the injected 
polaritons is higher (see fig. |D-F where no phase cut 
corresponding to a V-AV pair can be observed in the four 
regions delimited by the red squares). This confirms that 
polariton-polariton interactions, causing the widening of 
the subsonic regions and the consequent rearrangement 
of the polariton flows, induce the disappearance of V-AV 
pairs. 

IV. CONCLUSIONS 

We have investigated the washing out of vortex- 
antivortex lattices in exciton polariton systems as a func¬ 
tion of the polariton density. Our detailed analysis 
shows that V-AV annihilation and the washing out of 
the vortex lattice is due to polariton-polariton interac¬ 
tions that, through the renormalization of the lower po¬ 
lariton branch, induce the widening of the regions where 
the fluid is subsonic as the pump intensity is increased. 


This widening induce a modification in the polariton 
flows therefore changing the steady state position of vor¬ 
tices and antivortices until a threshold pump intensity is 
reached at which V-AV merging takes place. The com¬ 
plete washing out of the V-AV array can take place be¬ 
cause in our system has the number of vortices is equal 
to the number of antivortices, i.e. the system has zero 
angular momentum. Our analysis applied to a system 
with net angular momentum could open the way to the 
study of vortex-vortex interactions. 
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